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ABSTRACT 

We examine the distribution of injected energy in three-dimensional, adaptive-grid 
simulations of the heating of cooling flows. We show that less than 10 percent of the 
injected energy goes into bubbles. Consequently, the energy input from the nucleus is 
underestimated by a factor of order 6 when it is taken to be given by PVy/i^y — 1), 
where P and V arc the pressure and volume of the bubble, and 7 the ratio of principal 
specific heats. 
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1 INTRODUCTION 

It is now widely accepted that bipolar outflows from ac- 
tive galactic nuclei (AGN) largely offset the radiative losses 
of the hot plasma that forms "cooling flows" - systems in 
which there is a central depress ion in the temperature of the 
gravi tationaly trapped plasma (Bimicv 2005; Don ahue et al.l 
I2006J ). Notwithstanding this concensus about the impor- 
tance of AGN for cooling-flow evolution, considerable con- 
troversy surrounds the mechanisms by which AGN heat 

cooling flows. 

Observations, first by ROSAT feohringer et al.lll993l ) 



200lD of regions of 
enhanced syn- 



and later by Chandra (|Fabian et al 
depressed X-ray emission and (usually) 
chrotron emission, have played a crucial role in overturn- 
ing the paradigm of distributed mass dropout that domi- 
nated the field in the 1980s and 1990s. These regions have 
been c alled "bubbles" and considerable effort, bot h numer- 
ically (IChurazov et al.ll2000l: iBriiggen et ai~ll2002i) and the- 
oretic ally 1 Churazov et al. 20021 : iRovchowdhurv et al.ll20ol . 
I2005T) has been devoted to understanding the way bubbles 
interact with the surrounding, higher- density plasma. 

Since lines of sight through a bubble also pass through 
large volumes of X-ray emitting plasma in front of and be- 
hind a bubble, even a weak depression in the X-ray sur- 
face brightness at the location of the bubble may be indica- 
tive of a larg e drop in emissivity inside the bubble. In fact 
iFabian et al.l (|2000) report that the data for the Perseus 
cluster are consistent with zero emissivity inside one of the 
bubbles. Consequently, it is generally assumed that inside 
bubbles the density of thermal plasma is so low that the 
plasma emits negligibly in X-rays. This material is gener- 
ally assumed to be matter shot out from the nucleus in a 



E-mail: binney@thphys.ox.ac.uk; alouani@astro.ox.ac.uk 



high- velocity jet, which was thermalized to a high tempera- 
ture in the hotspot that marks the interface between the jet 
and shocked ambient plasma. 

Given that the plasma is close to an ideal gas, the ther- 
mal energy inside a bubble is U = PV/(7 — 1), where 7 is 
the ratio of principal specific heats. By taking the bubble to 
be in pressure balance with the X-ray emitting plasma im- 
mediately outside it, P can be determined from the X-ray 
data, while V follows from the geometry. Hence the energy 
inside bubbles can be reliably estimated from observational 
data. 

During inflation of a bubble work is inevitably done on 
the surroundings plasma. The minimum work that must be 
done at this stage is the work that would be done if the bub- 
ble expanded reversibly, which is PV, bringing the minimum 
energy associated with a bubble to 7/(7 — 1)PV = \PV for 
a non-relativistic plasma (7 = |). However, the inflation 
of bubbles is probably anything but reversible, and in this 
case, by Clausius' inequality the work done will exceed PV. 
We use high-quality hydrodynamical simulations of bubble 
formation to examine how much work is actually done by 
the jet, and conclude that it exceeds \PV by a substantial 
factor. 



2 THE SIMULATIONS 

The framework within which we si mulate cluster heat- 
ing has been describ ed elsewhere (|Omma et al.l |2004| ; 
lOmma fc Binnevll2004). We use the adaptiv e-grid Eulerian 
hydrocode ENZO (|Brvan fc Normanl Il997l ) with a piece- 
wise-parabolic Riemann solver. The plasma is assumed to 
be an ideal, monatomic, non-relativistic gas. Radiative losses 
are calculated in the optic ally-thin limit using the line and 
continuum cooling rates of ISutherland fc Dopital (|l993l ). 
The initial conditions are modelled on the current con- 
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Figure 1. The full curve shows the energy actually injected by 
the jets, while the broken curve shows the energy associated with 
bubbles. 



figuration of the Hydra c l uster as determined from Chan- 
dra data by iDavid et alJ <|200lJ >. Jets are launched by de- 
positing mass, m omentum and energ y in an injection disk 
of radius 3kpc l|Omma et alJ |2004| ). The grid cells are 
cubical and 0.6 kpc on a side at the centre, growing in 
stages by a factor of two in side length up to 36 kpc at 
the edge of the grid, which is a cube 614 kpc on a side. 
Periodic boundary conditions are imposed on this outer 
boundary. In contrast to many other simulations of cluster 
heating (e.g.. iBasson k, Alexander! l2003t iHeinz et al.ll200l 
IVernaleo fc Reynolds! [20061 ) there is no inner boundary. 

We hav e analyzed the energetic s of all the simulations 
described in lOmma fc Binnevl l|2004h . These differed in the 
amount of time the cluster cooled passively before jet igni- 
ton, and in the power and duration of the jets, as well as 
in whether there was more than one episode of jet ignition. 
However, we found similar energetics of bubbles in all simu- 
lations, so here we present only results from simulation S2, 
in which the cluster cooled passively for 300 Myr before a 
pair of jets with combined power 10 45 ergs _1 ignited and 
ran for 25 Myr. 



3 RESULTS 

In a series of snapshots bubbles were identified as regions 
in which the density was a factor / smaller than the mean 
density at that radius, and the volume V of the bubbles 
determined as a function of /. Clearly V is zero for / = and 
is comparable to the total available volume for f — 1. For 
0-2 ^ / ^ 0.75 there is a plateau in V that we identify with 
the bubbles, and we present results obtained by adopting 
/ = 0.3. 

In Fig. [T]the full curve shows as a function of time since 
jet ignition the energy injected into the simulation, while 
the dashed curve shows §PV for the bubbles. The energy 
associated with the bubbles peaks at 14 x 10 58 erg, while the 
energy injected is 78 x 10 58 erg. Thus less than a fifth of the 
injected energy is ever associated with identifiable bubbles, 
and at most epochs the fraction is less. 

To answer the question "where is the energy?" we have 
analysed the volume within 200 kpc of the cluster centre 
in the following way. We grouped the computational cells 
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ranges (0, r cavit y) (columns 2 to 5) and (r cavit y , r shock ) (columns 
6 to 9). The values of r cav j ty are 15, 37, 58 and 70 kpc, while those 
of r shock are 35, 75, 100 and 115 kpc. AE = AU + AK - AW. 



into approximately 40 spherical shells and in each shell i 
calculated 
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where V a is the volume of the cvth cell and e a is an energy 
density within that cell. The energy densities considered 
were u — ^nksT, where n is the particle density, w = — p$, 
where p is the matter density and $ < is the gravitational 
potential, and k = \pv 2 , where v is the flow speed. In this 
way we obtain measures AU, AW and AK of the changes 
since jet ignition in internal, potential and kinetic energy 
within each shell. 

Fig. [2] gives the radial distribution of energy differences 
at four times, 25, 50, 75 and 85 Myr after jet ignition. At 
25 Myr (top left panel) the jets have just switched off and the 
curves are structured by two distinct processes. At r < 30 kpc 
the jets are dominant, producing sharp peaks in internal (full 
curve), potential (short-dashed curve), and kinetic energy 
(dotted). Interior to the peak in potential energy AW at 
r ~ 20 kpc, there is a region in which AW < because 
uplift has reduced the mass density there. In this region the 
internal energy is only very slightly depressed because the 
drop in p is compensated by a rise in T. Beyond ~ 40 kpc the 
curves for AU and AW are structured by radiative cooling, 
which causes a slow decline in T and a flow of material 
from large to small radius that causes AW to be positive at 
r < 150 kpc and negative further out. 

The long-dashed curve in the top left panel shows the 
radiated energy R: unlike the other curves, this is not cal- 
culated by differencing energy densities at two epochs but 
by integrating the X-ray emissivity jx through the volume 
of each cell and with respect to time since jet ignition. 

The top right panel of Fig. [2] shows the disposition of 
energy 50 Myr after jet ignition. Compared to the top left 
panel for t = 25 Myr, the long-dashed curve of the radiative 
losses has risen at all radii, and the peaks in AU, AW and 
AK have moved outwards with the rising jet debris. Now the 
curves for AU, AW and AK are structured by the jet for 
r < 70 kpc and by cooling further out. At small radii the peak 
in AU is now preceeded by a pronounced trough, reflecting 
the development of an under-pressured region behind the 
zone filled by material that is being carried up by its mo- 
mentum, and decelerating in the cluster's gravitational field. 
At r > 80 kpc cooling continues to drive an inward flow that 
makes the negativity of AU and AW at large radii more 
pronounced. 

The bottom left panel of Fig. for t = 75 Myr after 
jet ignition, shows a new feature: the curves for AU and 
especially AW have developed peaks at r ~ 20 kpc, interior 
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Figure 2. The radial distribution of energy changes defined by equation JTJ since jet ignition. 
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to the troughs that preceed the main peaks. This new region 
of enhanced potential energy reflects material that has fallen 
down after being uplifted and is now at a radius lower than 
the radius it would have in an equilibrium model, given its 
specific entropy. Thus this peak reflects a particular phase 
of the gravity waves that are strongly excited by the jets. 
The bottom-right panel, for r = 85 Myr is very similar to 
that for 75 Myr. 



At each panel of Fig.[2]one can identify two characteris- 
tic radii, namely the radius r s h oc k that marks the end of the 
sharp drop in curves for AU and AW that the shock pro- 
duces, and the radius r cav ity at which these curves cross zero 
on the other side of the peak. For the time of each panel in 
Fig. [2] Table[T]gives energies, in units of 10 58 erg, contained 
in the regions r < r cav ity and r- cav it y < r < r shock . We see 
that at 25 Myr, as the jet shuts off, the total energy inside 
^cavity has increased by an amount AE — AU +AK—AW = 
123 x 10 58 erg, that greatly exceeds the total injected energy, 
£j c t = 78 x 10 58 erg and is more than an order of magni- 
tude greater than the energy then associated with bubbles 
(8.9 x 10 58 erg from Fig. [TJ . This energy increase is largely 
caused by a negative value of AW = —132 x 10 58 erg fol- 
lowing uplift of material - to achieve this value, the jet has 
not in fact supplied 132 x 10 58 erg of energy, which would 
suffice to lift the uplifted material to r — oo, but only the 
smaller amount required to move material out of the cav- 
ity zone. This explains why AE for this region exceeds the 



the jet energy. In the dense region outside r. 



cavity 



and inside 



r s hock, AW has increased by twice the jet energy, reflecting 
(i) the arrival from below of gravitationally bound uplifted 
material, and (ii) the arrival from above of material that 
has cooled. The internal energy of the gas in this region has 



risen by 75% of -Ejet- At this stage the bulk kinetic energy 
alone is equal to the energy associated with bubbles. 

At later times we see that within r cav it y , AU and AK 
continue to increase, while AW becomes less negative. In the 
dense shell further out AK continually decreases, while AU 
slowly rises. The total energy change AE in the dense shell 
is always falling as material is driven in by the jet-induced 
flow on one side and gravity on the other, bringing with 
it large (negative) potential energies. In this region cooling 
also attenuates the rise in AU that would otherwise occur. 



4 DISCUSSION 

We have shown that in the simulations the minimum energy 
associated with bubbles, _E m i n = §PV, is smaller than the 
injected energy by at least a factor of nearly 6. Why is this 
so, and can we expect the result to extend to real cooling 
flows, or does it merely reflect limitations of the simulations? 

The factor | allows for the energy %PV that is actu- 
ally in the bubble, and an allowance PV for the work that 
must be done to create a cavity in the ambient medium. 
As explained in the Introduction, this is the minimum al- 
lowance, and is the actual requirement only if the cavity is 
created reversibly and at constant pressure. In this case the 
displaced material would be neither heated nor accelerated, 
and the work PV would be done on the mechanism (piston 
etc) that held the pressure constant. In reality the creation 
of the cavity is strongly irreversible, so the displaced mate- 
rial is heated and accelerated. In as much as real jets prob- 
ably involve larger velocities than those (~ 40000 kms -1 ) 
characteristic of the simulations, the simulations may be ex- 
pected to under- rather than over-estimate the irreversibility 
of bubble inflation. 



4 J. Binney, F. Alouani Bibi & H. Omnia 



The internal energy of displaced material is changed 
by of order itself, and its original value was equal to the 
bubble energy. Moreover the material evacuated from the 
cavity itself displaces material, so the mass of the material 
that has to be moved and heated exceeds the mass originally 
associated with the cavity by a factor a > 1, say. Doubling 
the internal energy of the shifted material, would involve 
supplying |<xPV of energy over and above what is required 
in the reversible case. 

In the simulations the timescale for bubble inflation is of 
order the sound-crossing time at the ambient sound speed, 
and it is widely assumed that the same is true in real clus- 
ters. Consequently, the shifted material inevitably gains an 
amount of kinetic energy that is comparable to the original 
internal energy, ^aPV, as Table [1] confirms to be the case. 

These two estimates of energy that must be supplied 
should be added to the energy "/PV/(-y — 1) required in the 
reversible case, where 7 is the ratio of principal specific heats 
of the material in the cavity. In our simulations 7 = | , so the 
actual energy required is larger than the minimum energy 
by a factor 1 + |tv. If the cavity contains relativistic plasma, 
as is frequently assumed, and the ambient plasma is non- 
relativistic, the energy that must be supplied is larger than 
the minimum energy by a factor 1 + |q. The energetics of 
the simulations implies that the shifted mass is a = 3.75 
times the mass originally in the bubble 

In reality a is probably not so large because diffusivity 
will lead to PV being underestimated, both in the simula- 
tions and in real cooling flows. Our simulations do not ex- 
plicitly include thermal conductivity, but any numerical hy- 
drodynamic scheme has significant diffusivity on the scale 
of the grid (or the mean-free path in the case of smooth- 
particle hydrodynamics) by virtue of assigning single val- 
ues to the dynamical variables at each grid point, when in 
reality the material in a cell would contain hot and cold 
regions. As a consequence, energy can cross a cell in one 
timestep when hot matter enters on one side at the start 
of a timestep, and warmed matter leaves on the opposite 
side at the start of the next timestep. There is no concen- 
sus regarding the appropriate value s of the kinetic tr a nspor t 
coefficients of int racluster gas [see iKim fc Naravanl l|2003h ; 
iPope et al.l (|2005T ) for conflicting positions], and even less un- 
derstanding of the effectiveness of turbulent transport. How- 
ever, turbulence stretches and distorts bubbles, and when 
their narrowest dimension becomes comparable to the grid 
spacing, diffusivity will eliminate them. Thus kinetic and 
turbulent diffusivity cause bubbles to shrink on a dynamical 
timescale in the simulations, as they probably do in reality 

The broken curve in Fig. [1] illustrates this phenomenon 
at two points: between t ~ 8 and 14Myr, PV dips while the 
jet is still firing, and then at t > 50 Myr it falls continuously 
after the jet has switched off. The dip at early times reflects 
shrinkage in the volume of low-density material by ~ 25 per- 
cent following the emergence of turbulent eddies around the 
jet; as the jet establishes itself at t m 7 Myr, the cavities shift 
from being nearly hemispherical to elongated (Fig. [3j . The 
steep velocity gradient perpendicular to the jet axis then in- 
duces large-scale turbulence, which draws cool material into 
the cavity. The ensuing fragmentation of the cavity leads 
to the decline in cavity volume that causes the dip in the 
PV curve of Fig. [T] At t > 14 Myr a new regime is entered 
in which lobes swell around the now more slowly advancing 




Figure 3. Cross sections through the cavities at eight times. In 
addition to the time since jet ignition, each panel gives the total 
volume of the cavities. Note that there are scale changes between 
images. 

hotspots. The drop in the PV curve at t > 50 Myr is also 
driven by turbulent breakup. How realistic the steepness of 
this drop is, we do not know. 

Observationally it is inevitable that the volume of cav- 
ities will be underestimated because it has been determined 
from X-ray surface brightness, rather than from a knowledge 
of the plasma density at each point in three dimensions, as 
can be done in the simulations. Clearly such underestimation 
of the volume will lead to a corresponding underestimation 
of the power input by the jets. It would be interesting to 
compare estimates of cavity volumes obtained from X-rays 
with ones obtained from high-resolution radio maps, which 
might be less sensitive to projection effects. On the other 
hand, the X-ray cavities sometimes do not coincide with the 
radio lobes - the Hydra cluster p rovides an example of this 
phenomenon l|Nulsen et al.ll2005T ). 

This discussion leads to the conclusion that taking the 
energy injected by jets to be \PV probably does lead to its 
being underestimated by a substantial factor, which we have 
estimated to be ~ 5.5. Our estimate of this factor may be on 
the high side if the simulations over-estimate the diffusivity 
of the plasma that arises from both turbulent mixing and 
kinetic processes, but may be on the low side because our 
jet speeds are relatively low and cavity volumes are likely to 
be underestimated from X-ray surface brightnesses. 

The fact that most of the injected energy is never in 
cavities, has important implications for the time- and space- 
dependence of cluster heating. As Table [1] shows, much of 
this energy is stored as potential energy - the energy re- 
quired to lift low-entropy gas above its natural position 
in the cluster atmosphere. This energy subsequently heats 
the plasma through the dissipation of the non-linear grav- 
ity waves to which it gives rise. The timescale of this re- 
lease can be very long, especially in the case of an unusu- 
ally violent outburst of the AGN, which drives cool plasma 
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to great heights, at which the dynamical time can be long 
jOmma et al.ll2004 ). 

Several s tudies 
Nipoti fc Binnevl 



120051 : 



(e.g., 


3irzan ct 


all 22041: 


iBest et al. 


20061; 


lAllen et al.1 



20061 ) have considered the implications of cooling flows for 



the phenomenolgy of AGN and radio sources under the 
assumption that the injected energy is yPV/^ — 1). The 
studies conclude that heating by AGN can balance radiative 
cooling provided the duty cycle of the jets is high. With our 
increased estimate of the energy injected in a given episode 
of bubble inflation, the required duty cycles will be shorter. 

Assuming a 10 perc ent efficienc y in th e conversion of 
rest mass to jet power, lAllen et al.l |2006) conclude that 
jet powers are smaller than the rate of energy release from 
Bondi accretion of hot plasma by a factor ~ 6 (their Fig. 4). 
Prima facie, our increase in the required jet power would 
bring it up to the rate of release of energy by Bondi ac- 
cretion. Is this conclusion problematic, given that only a 
fraction of the plasma that passes over the Bondi accretion 
radius can be swallowed by the black hole? 

The rate at which material passes the Bondi radius pro- 
vides a useful measure of the mean jet power, but not of the 
instantaneous power, because inside the Bondi radius an ac- 
cretion disc is likely to provide a reservoir of material, and jet 
outbursts ar e thought to be b r ief episodes in the ev olution 
of this disc (jNipoti et al.ll2005l ; iKoerding et al.ll2006r ). Thus 
on the timescale of bubble inflation, jet power can safely ex- 
ceed the mean power available by Bondi accretion, even if 
the black hole's only mass supply is from virial-temperature 
gas (which is by no means certain) . Estimates of the power 
available from Bondi accretion are most use fully compare d 
with the radiative losses of the cooling flow l|Binney||2005f ). 
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